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Abstract. We consider Dirichlct eigenfunctions ux of the Bunimovich sta- 
dium S, satisfying (A — X^)ux = 0. Write S = RUW where R is the central 
rectangle and W denotes the "wings," i.e. the two semicircular regions. It is 
a topic of current interest in quantum theory to know whether eigenfunctions 
can concentrate in i? as A — ► oo. We obtain a lower bound CX~^ on the 
mass of ux in W, assuming that ux itself is L^-normalized; in other words, 
the norm of Ux is controlled by times the norm in W. Moreover, if 
Ux is a o(A~^) quasimode, the same result holds, while for a o(l) quasimode 
we prove that norm of Ux is controlled by A* times the norm in W. 
We also show that the norm of ux may be controlled by the integral of 
to|9jviip along dSn W , where w \s a. smooth factor on W vanishing at Rn W . 
These results complement recent work of Burq-Zworski which shows that the 
norm of ux is controlled by the norm in any pair of strips contained in 
R, but adjacent to W . 



1. Introduction 

The Bunimovich stadium S* is a planar domain given by the union of a rectangle 
R — {{x,y) I X G [—a, a], y G [—/?,/?]} with two "wings," i.e. the two semicircular 
regions centered at (±a, 0) with radius /3 which lie outside R. Geodesic flow in 
S (obeying the law of reflection at the boundary) was proved to be ergodic by 
Bunimovich pQ. It follows from this and from results of Schnirelman 0, Zelditch 
and Colin de Verdiere [3 that the stadium is quantum ergodic. This means 
that there is a density one sequence of Dirichlet eigenfunctions which becomes 
uniformly distributed; in particular, along this density one sequence the weak limit 
of the mass distribution becomes uniform. One can ask whether the entire 
sequence of eigenfunctions becomes uniformly distributed; if so, the domain is called 
quantum unique ergodic (QUE). It has been conjectured by Rudnick and Sarnak 
that complete surfaces with negative curvature which are classically ergodic are 
QUE; this has been proved recently by Lindenstrauss 6 for arithmetic surfaces.^ 
The Bunimovich stadium, by contrast, is generally believed to be non-QUE; it is 
thought that there is a sequence of eigenfunctions that concentrates in the rectangle 
R. Little is currently understood about the way in which such eigenfunctions would 
concentrate, however. For example their (hypothetical) rate of decay outside R is 
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unclear. The result in the present paper is intended to shed some light on this 
question: we show that any sequence of eigenfunctions (or quasimodes) cannot 
concentrate very rapidly inside R, by obtaining lower bounds (tending to zero as 
A — > oo, but only polynomially) on the mass inside the wings W±. 

Let A = —d^ — dy denote the (nonnegative) Laplacian on S with Dirichlet 
boundary conditions. We denote by ||-|| the norm in L'^{S), and by dpfg the outward 
pointing normal derivative of g at dS. We consider a o(l) Dirichlet quasimode ux 
for A, by which we mean that we have a sequence A = Afe ^ cxd of real numbers 
and a corresponding sequence u\ e H'^{S) satisfying 

(A - X^)ux - /a, 

(1) ux las = 0, 

II«a|| = 1, 

where 

(2) |l/A||=o(l)asA^^. 

We more generally define a 0(A^^) or o(A^^) quasimode by modifying the right- 
hand side of (0) accordingly. Of course a sequence of eigenfunctions is a o(A~^ ) 
quasimode for any j. 

It is easy to see that a 0(1) quasimode can be localized to a small rectangle of 
the form [7,(5] x [— /3,/3], where [7,(5] is an arbitrary subinterval of [— a,Q;]; indeed 
the family U(„+i/2)7r/,9 — (f>{x)cos{{n + l/2)7rj///3) is (after normalization) such 
a quasimode, where (p is any nonzero smooth function supported in [7,(5]. By 
constrast, an o(l) quasimode cannot be so localized: Burq-Zworski have shown 
that the norm of ux is controlled by (that is, bounded above by a constant times) 
its norm in the union of any two rectangles of the form ([—a, 71] x [— /3,/3]) U 
([72, a] X [— /3,/3]). In particular, for a o(l) quasimode, the mass cannot shrink 
to a closed region disjoint from the wings of the stadium as A ^ 00. 

Although the stadium is classically ergodic, there is a codimension one invariant 
set for the classical flow, consisting of vertical "bouncing ball" orbits parallel to the 
y-axis and within the rectangle R, and the union of these orbits is the most likely 
place where localization of eigenfunctions, or more generally o{l) quasimodes, can 
occur^. There is a rather convincing plausibility argument in the physics literature 
due to Heller and O'Connor which indicates that a density-zero sequence of 
eigenfunctions, with eigenvalues ((n -I- l/2)7r//3)^ -I- 0(1), does concentrate to some 
extent at these bouncing-ball orbits. The rigorous essence of this argument has 
been developed by Donnelly ^ who showed that there are sequences of functions 
lying in the range of spectral projectors Ej^{A), where /„ are intervals of the form 
[{{n + l/2)7r//3)2 - C, ((n -I- l/2)n/f3)^ + C] which concentrate at the bouncing- 
ball orbits.^ On the other hand, the result of Burq-Zworski [2j shows that such 
localization cannot be too extreme: the control region must extend to the boundary 
of the rectangle. 

Our main result here is that we may in fact push our control region outside the 
rectangle altogether and into the wings, in return for a loss, either from restriction 
to the boundary, or of powers of A. To state this concisely, it is convenient to 

^The explicit quasimode in the paragraph above concentrates along a subset of these orbits 
■^This was shown for surfaces without boundary containing a flat cylinder, but the arguments 
go through for the stadium. 
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introduce an auxiliary coordinate in the wings given hy w — \x\ — a; thus w is 
nonnegative on the wings and vanishes exactly on the vertical lines RDW. 

Theorem 1.1. There is aC > 0, depending only on a/ (3, such that any family ux 
satisfying obeys the estimates 

(3) ||/a||'+ / w\dNUxfdl>C , 

Jdsnw 

(4) ll/Af + A«||^A|li.(vK)>C 

and 

(5) \'\\fxr + XHnx\\l.^^^>C. 

Therefore, if ux is a o(l) quasimode, we have for sufficiently large A 

\\w^dNu\\L'^(^gSnW) > C, 

(6) 

\\ux\\l2(w)>(^^ ' 
while if Ux is a o(A^^) quasimode (e.g. an eigenf unction), 

(7) hxWmw) > c\-^ 

Note that the results of Theorem ll.ll still leave open the possibility of quasimodes 
concentrated along bouncing-ball orbits in the rectangle with o(l) mass in the 
wings. They also do not rule out the possibility that all the energy in the wings 
may asymptotically concentrate in a boundary layer near R. 

2. Preliminaries to estimates 

Our main tool is positive commutator estimates, which we use in the following 
form: 

Lemma 2.1. Let u he real, equal to zero at dS, and satisfy (A — X^)u — f, where 
f is smooth. Then for any real vector field A, 

(8) {u,[A- \^,A]u} = ({2Au+ {divA)uJ} + [ {dNu)Audl. 

JdS 

Proof. We integrate twice by parts, using the Dirichlet boundary conditions in the 
first instance, to write 

{u,[A- X^,A]u) = {f,Au) + I dNuAudl- {u,Af). 

JdS 

Applying Green's Theorem to the last term now gives two terms: {Au, /)-|-((div A)u, /). 
Since u and / are real this yields the desired identity. □ 

We also record here an inequality that will be of use in estimating derivative 
terms. 

Lemma 2.2. Let u, f be as in Lemma \2.1\ 

Then for all s > 0, for X sufficiently large, 

\\Vuf < Cs{\"'^''^^-'^\uf + \-'\\ff). 
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iVwII^ = / ul+uldA 



Proof. We compute 



{Au)udA 

[ u^dA+{f,u). 
Js 

Applying Cauchy-Schwarz to (/, u) gives the estimate. □ 

3. Proof of © 

It suffices to prove ^ under the assumption that ux, and hence /a, are real, since 
we can treat the real and imaginary parts separately. We make this assumption 
from now on. 

We begin with the standard commutator [A,a;9a:] = —29^. Applying Q with 
A = xdx, and dropping the subscript on u\, we have 

(9) f f 

= / xdxudNudl + / {2xdxU + u) f dA; 

JdS J 

in the last equation we integrated twice by parts, using for a second time the fact 
that (A — A^)m = / as well as the fact that u satisfies Dirichlet boundary conditions, 
hence integration by parts produces boundary terms only where derivatives land 
on both factors of u. 

Now at every boundary point we may decompose xdx into pdi + qOn where di is 
differentiation tangent to the boundary. Of course di annihilates u. Now since dx is 
tangent to the upper and bottom sides of the rectangle we find that the boundary 
integral in ^ is only over dSnW. Moreover, as is tangent to the top and bottom 
of the circles forming the boundaries of the wings, we have q = 0{w) on dS fl W. 
Hence we have shown that 

(10) ||u:.f</ 0{w)\dNufdl + e [{u^ +ul)dA + C [ fdA. 

JdSnW J J 

We may absorb the e||ua;|p term, then apply the Poincare inequality, and absorb 
the e||u|p to obtain 

(11) \\u\\l.,s.<C [ w\dNufdl + C\\ff 

JdSnw 

which is the first part of our theorem, as we took u to be L^-normalized. 



4. Proof of Q 
To prove this estimate we start from 



(12) \\u,r < c / w+\dNu\'di+c\\fr 

JdS 

which follows directly from the considerations of the previous section, and estimate 
the boundary integral term. We shall obtain upper bounds of the form 



/ u^dA+WfW 
Jw 



A 

iw 



SPREADING OF QUASIMODES IN THE BUNIMOVICH STADIUM 



5 




Figure 1 . The three regions of interest in W. 



and 

/ u^dA + X'WfW^, 
Jw 

thus proving Q and 

We shall perform this estimate in three separate regions in the wing. Region I 
is the near-rectangular region, in a boundary layer where w < SX~^. Region II will 
be outside the boundary layer, where 5X~^ < w < f3/2. Region III will be the far 
outer region w > (3/2. 

We begin with Region III, far away from the rectangle. In this case we employ 
Lemma |2 . 1 1 where A is the operator (f)(x)dj:, where is supported where w > /3/4, 
say, equal to 1 where w > (3/2, and with dx(t> > 0. Then (jHJ gives us, with 
P^[^,A], 

/ (jjd^udMU < \{Pu,u) \ + |(0x/,it)| + 2{(j)UxJ). 

JdS 

Note that P = -24)xdl^ - (p^A- Thus the LHS is bounded (Ve > 0) by 

\{-(j)xUxx,u)\ +<i\\Ux\f +C{\\u\\\2^w) + ll/f ) 

(where of course C depends on e). We can add the positive term Jg <px{dyu)'^ to 
this estimate. Integrating by parts in y gives us 

\{M-Uxx-Uyy),u)\ +e|lU:rf +C(||w||^2(vv) + ). 

Using the positivity of the integrand, we thus obtain an estimate 

\2|| ||2 , rn2 , ,, ,,2 



(f> dxU OnU dl 

asniii 



<CX'\\u\\l.^^^+C\\fr + e\\uxr 



(13) 

with C depending on e > 0. 

Now we work on Region I, within a 0{X~^) boundary layer along the rectangle. 
We again apply Lemma [2. II this time with A = xdx + ydy. Since A is a tangential 
vector plus a positive multiple of On all along dS we obtain 

{dNufdl < \{u, [A- X'^,A]u)\+2\{uJ)\ + 2\{AuJ)\. 

as 

Using [A — A^, ^] = 2 A and Cauchy-Schwarz this becomes 

\dNu\^dl<CX^\\u\\^ + \\Auf + C\\f\\^. 

dS 
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Restricting to 35 fl I in the integrand, we can estimate w in by (5 A ^, and this 
gives 

(14) / w+ldNufdlKSCiWuf + X-^iWff + WAuf)). 

Jdsni 

Using Lemma f2. 21 we may estimate by C(A^||w||^ + ||/||^). Hence we may 

finally write 



(15) / w+l^jv^rd/ <(5Co||u|r + (5CA-^||/|r; 

Jdsni 

note that in the above construction, Co can in fact be chosen independent of S. 

Finally, we estimate in Region II. To begin with, we note that for w > 6X~^, we 
can estimate w+ by (5~^A^w^, so we have 

/ w+\dNuf dl < CS^^ / }?w\x{v)dyudNudl\ 
JdSnii JdSnii 

here we take x supported in \y\ > (3/20, and equal to —1 for y < —(3/10 and +1 for 
y > (3/10, so that x{y)(^y is a positive multiple of plus a tangential component 
on dS n II. To estimate further, we employ Lemma l2. II with the commutant 



A = x'wixiy)d. 



V 



The point of this commutant is that we have given ourselves two powers of w which 
will "absorb" two integrations by parts in x without any boundary terms at w = 0; 
this is crucial since we know of no way to deal with such boundary terms. (On the 
other hand, we pay the price of additional powers of A with this gambit.) Thus we 
obtain, setting Q = [A, 



(16) / X^wXx{y)dyudNudA< {Qu,u) 
Jas 



+ X' 



w+xyiy)uf 

s 



2|(A«,/}| 



We can estimate the second term on the RHS by A'*||u||^2(^;i/) + C||/p. 
Now consider the terms involving the operator Q. This is given by 

Q = A2 ( - 4w+xd^dy - 2wlxydly ~ 2H{w)xdy - wlxyydy) 

where H{-) is the Heaviside step-function, H{w) = for w < and 1 for w > 0. 
To treat the terms involving one derivative, e.g. the third term above, we integrate 
by parts: 

-2A2 / H{w)x{y)dyUu = X'' f H{w)xy{y)u^ 

which is therefore estimated by A^ ||u|||2(vi^-) • The fourth term is estimated in exactly 
the same way. 

Thus we are left to estimate 

(17) A^ ( - 4{w+x{y)dlyU, u) - 2{w\xy{y)dlyU, u)) 

Integrating the first term by parts in x and the second term by parts in y gives us 
two principal terms 

(18) 4A^ (^{wll_'^dxU, wll'^dyu) + {w\xydyU, dyu)^ 
together with two other terms 

AX^(^{H{'w)dyU,u) + {w\xyy{y)c)yU,u) ^ 
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which are estimated in the same way as the first order terms above. We apply 
Cauchy-Schwarz to the first term in H18|l . while in the second term, which is positive, 
we replace by (which is larger, up to a constant multiple) and then integrate 
by parts again, getting (up to another first-order term estimated as above) an upper 
bound for H18|) of the form 



(19) CX' / w+\d^u\" + w+\dyu\" dA. 

Js 

Now we integrate by parts again, getting 

(20) CA^y (^w+{—Uxx — 'U'yy)u + H{w)uUx^dA. 
Writing ~Uxx — Uyy = X?u + /, we can estimate the integrand of H2Uf) by 

A'*u'+|m|'^ + X^w+fu + \^ H{w)uux 

< X'^w+\u\^ + ^{X^w+u^ + w+f) + CX'^H{w)u^ + eul- 

Hence we may estimate (|18f) . by as 

CX''\\u\\l.^^^+e\\ux\f + C\\ff. 

The term 2\(Au, f)\ is bounded in a similar manner: By Cauchy-Schwarz we 
may estimate it by 



or by 



C\\ff + CX^w+Uyf 



CX^\\ff + CX'\\w+Uyf, 



as we prefer. Our estimate for l|19() turns the latter estimate into 

CX\ff + CX^u\\l,^^^+e\\uxf. 

On the other hand, treating A** || w+UyH^ in the same manner gives a bound by 

C\\ff + CX'\\u\\l^^^+e\\uxf. 

On the support of x ^-iid at the boundary, x{y)^y is a positive multiple of d^u. 
So the upshot is that (|16|l now gives 

(21) / X'wlx{y){dNuf <CX'\\u\\l,^^^+e\\uxf + C\\ff. 

JdS 

and 

(22) / X^wlx{y){dNuf<CX^u\\l,^^^+e\\uxf + CX^\\ff. 

J dS 

At last we can estimate the boundary integral term in (|12|l by a combination of 
lUSJ), (O, and (EU/l|221), obtaining 

w+\dNu\^ dl < SCoWuf + CX^\\u\\l,,^. + 2e\\uxf + C\\ff. 

as 
and 

/ w+\dNu\' dl < SCoWuf + CA^Il^ll^ + 2e\\uxf + CX'\\ff. 

JdS 
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Here C depends on e, 5 but Co is independent of 6. We now combine this with l|12|l . 
Absorbing the \\ux\\ and ||m|| terms on the LHS by taking S and e sufficiently 
smah, we obtain (O. 

5. Concluding remarks 

The estimates presented here are certainly not optimal. The powers of A appear- 
ing in Theorem ll.f l ean probably be improved using refinements of the methods used 
here, but it seems unlikely that one could achieve an optimal result with them. We 
have not, therefore, attempted to obtain the best possible powers of A, but have 
rather tried to present a poynomial lower bound on |jwA||L2(^^) with a simple proof. 

It would be of great interest to obtain a polynomial lower bound on the mass 
of u\ in a subregion of W which is a positive distance from R (i.e. region III of the 
previous section). We do not know whether such a bound holds, but it does not 
seem to be obtainable using the methods of this paper; possibly it might yield to 
the use of more sophisticated tools from microlocal analysis. 
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